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ABSTRACT 

We  consider  an  epidemic  model  of  the  form  S  -*■  I  S  with 
history  on  (-°°,0].  The  well-known  threshold  phenomenon  is 
discussed  in  terms  of  the  stability  of  a  functional  difference 
equation,  also  known  as  the  translation- invariant  renewal  equation. 
Since  the  difference  equation  has  infinite  delay,  the  work  of 
other  authors  on  finite-delay  problems  is  extended.  Also,  epidemic 
models  with  spatial  effects  are  discussed  by  extension  of  the 
results  to  difference  equations  in  a  Banach  space. 
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FUNCTIONAL  DIFFERENCE  EQUATIONS  AND  AN  EPIDEMIC  MODEL 


Recently  there  has  been  interest  in  a  model  for  the  evolution 


of  a  disease  given  by  the  equation 


SC1)  "  S(t)  |°  B(0)  ^  S(t+8)d0, 


where  S(t)  is  the  density  of  susceptible  individuals  and  the 
right  hand  side  allows  for  infectious  contacts  at  the  present  time 
due  to  the  past  history  of  infections.  Equation  (1)  is  an 
S  ■*  I  model,  that  is,  it  allows  only  for  susceptibles  to  become 
infected.  Diekmann  [2,3]  has  considered  this  model  and  also  has 
allowed  spatial  effects.  In  an  S  -*■  I  S  model,  where  one  allows 
for  recovery  (without  immunity)  of  infecteds,  one  would  add  a  term 

f°  i 

j  c(0)  S(t+9)d9  to  the  right-hand  side  of  equation  (1). 

J  —  oo 

The  "threshold  phenomenon"  of  Kermack  and  McKendrick  is  well- 
known  in  mathematical  epidemiology;  see,  for  example,  Bailey  [1] 
or  Hoppensteadt  [11].  For  equation  (1)  the  appropriate  initial  data 
consists  of  the  value  k  =  Sq  =  S(0)  and  the  initial  history 

<l>(0)  51  '35'  S(0),  9  €  (-<*>,0].  The  epidemiological  model  is  realistic 
only  if  SQ  >  0,  <j>(9)  <  0  and  B(9)  >  0  for  9  £  (-°°,0],  and 

B (.•)$(.')  t  0.  With,  these  assumptions  there  always  exists 

,Sq)  de=^n  s«oCt;SQ).  The  threshold' phenomenon  can 
be  stated  as:  There  exists  an  S*  such,  that 


(a)  For  fixed  SQ  >  S* ,  lim^S^C*  ;SQ)  >  SQ. 

(b)  For  fixed  SQ  <  S* ,  lim^S^C*  ;SQ)  -  SQ. 


2 


Here,  lining  stands  for  the  limit  in  the  function  space  where 
the  initial  history  comes  from. 


The  equation  (1)  can  be  transformed  into  an  equivalent 

n  for  y(t)  d®fn  ^  S(t)  by  nc 
y(T)dt.  Specifically,  (1)  is  equivalent  to 


functional  difference  equation  for  y(t)  ^  S(t)  by  noting 


that  S (t)  =  Sq  .+  | 


0 


(2) 


y  (t)  -  k  |  B(0)y(t+9)d0  =  Q  y(t+0)d0^Q  B(0  )y  (t+0  )  d0^ , 


where  we  have  notated  k  =  Sq  so  as  to  distinguish  Sg  as  a 
parameter  to  be  adjusted.  We  will  show  that  the  zero  function 
<(>(0)  =  0  for  0  €  (-°°,0]  is,  for  equation  (2), 

(a')  uniformly  asymptotically  stable  for  k  <  S* 

(b ' )  unstable  for  k  >  S* 

(C)  s*  =  i/(|°  B(0)de), 


and  then  we  will  show  how  the  stability  results  (a'),  (b ' )  for 
equation  (2)  imply  the  threshold  phenomenon  (a),  (b)  for  equation  (1) . 
Thus,  we  will  show  that  the  threshold  phenomenon  is  a  particular 
case  of  stability  results  for  functional  difference  equations. 

The  greater  part  of  the  mathematical  analysis  of  this  paper 
will  be  the  discussion  of  a  particular  class  of  functional 
difference  equations  which  will  be  wide  enough  to  obtain  the  result 
of  the  threshold  pehnomenon.  We  will  consider  the  equation 


9 


where  y  is  a  function  of  9  £  (-",0]  defined  by  yt(9)  =  y(t+9) 
and  D  and  G(t,0  are  defined  on  function  :  (-",0]  -+lRn,  with 
r  0  yQ 

D<j>  =  (0 )  -  e  A(9)<|)(9)d9  so  that  D  is  linear  and  of  a  re- 

J  _  oo 

stricted  nature.  Here,  A(9)  is  lRn><n  (real  n*n  matrix) -valued . 

The  positive  real  number  Y  is  assumed  to  be  fixed  and 

r0  n 
|  |A(9)|d9  <  ",  where  |*|  is  the  Euclidean  norm  for  R  and 

J  -oo 

the  corresponding  matrix  norm  on  lRn  n.  By  consideration  of  the 
spectral  theory  developed  below  it  will  hopefully  become  clear  why  the 

Y0 

special  form  e  A(0),  with  ) A (9 ) | d©  <  ",  is  assumed  for  the 

J  -OO 

kernel  of  the  integral  term  in  the  linear  operator  D. 

Since  (3)  has  infinite  delay  some  care  must  be  taken  in  the 
assumptions  made  about  the  class  of  equations  and  class  of  solutions 
to  be  considered.  Infinite  delay  in  retarded  functional  differential 
equations  has  been  considered  by  several  authors,  including  Hale  [5], 
Hale  and  Kato  [7],  and  Naito  [12].  We  will  look  for  solutions  yt 
in  the  space 


UCI11  -yfl 

CY  =  {4>:  (-",0]  -*-]R  |  $  is  continuous  and  there  exists  lin^^e 


With  D  restricted  to  the  above  special  form  it  would  not  be 
difficult  to  work  in  the  space  Kn  *  L  Y ,  where 

P  f  » 

,0 

®n  x  Lp>y,  where  Lp>Y  -  {*|j  ePyQ  (.9)  |Pd9  <  «}, 


but  we  will  not  discuss  this  further.  Also,  it  should  be  possible  to 


r 
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develop  a  theory  for  the  abstract  spaces  of  Hale-Kato  type,  as  in 
Hale  and  Kato  [7]  and  Naito  [12]. 

Note  that  for  every  4>  €  we  have  (a)  00  >  ($1^ 

^e=n  sup0 <0eY9  I ^  (e )  I  >  and  (b)  the  function  eY‘<(>(-)  is  uniformly 

continuous  for  9  £  (-",0].  is  a  Banach  space  when  given  the 

norm  | ■ |  .  Because 

I  D4>  |  <  (l  +  J°  |A(0)|d0)  |<f>  |  Y, 

D:  ]Rn  is  a  bounded  linear  operator. 

For  the  general  initial  value  problem  (4)  with  initial  data 
yQ  =  <P  €  we  allow  any  D:  -►  ]Rn  of  the  form 

f°  y8  r° 

D <p  =  <P  (0)  -  er  A(0)<f»(0)d0  with  |A(9)|d0  <  °° 

J  -oo  J  -00 

and  we  allow  any  G  €  C1QR+  x  C^;  lRn)  for  which  the  Frechet  differential 
differential  D^G(*,*)  is  bounded  on  all  sets  1R+  *  B,  B  any 
bounded  subset  of  C^.  For  any  initial  data  <P  €  Cy  fox  which 
EW  =  G CO ,4> ^  there  exists  a  unique  global  solution  yt,  i.e., 

Dyt  =  G(t,yt)  and  limt^Q  jyt-<$|y  =  0,  as  can  be  proven  using  the 
contraction  mapping  theorem. 

When  G  2  0,  so  that  (4)  is  linear,  there  is  a  unique  global 

solution  yt  *  yt(4>)  for  all  initial  data  <j>  satisfying  D<j>  =  0. 

Defining  TD(t)$  *  yt»  we  find  that  Tp(t)  is  a  semi-group  of 

defn 

bounded  linear  operators  on  the  Banach  space  SB  »  e  C Y :  Dd»  =*  0 > 

with  the  |.|  norm.  One  can  show  that  TDCt)  is  generated  by 
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the  operator  Ap  defined  by  AD<t>  =  <P ,  *  ~  gg"  *  on  its  domain 
^(Ad)  =  {$  £  <P  £  2). 

The  spectral  properties  of  Tn(t)  can  be  found  by  examining 

"  0 

the  characteristic  equation  0  =  det  A(X)  =  detf I  -  f  e^  +  Y^  A(9)d0 

'  *  -00 

as  is  made  rigorous  by  considering  the  generator  A^  and  the 
particular  nature  of  the  equation  Dyt  =  0.  Following  Naito  [12], 
let  us  notate  I  *  {X:  Re  X  >  -Y}.  Given  a  linear  operator  B, 
we  say  that  X  is  in  (a)  the  resolvent  of  B  if  (XI-B)  *  is 
well-defined  and  bounded, (b)  the  spectrum  of  B  if  it  is  not  in 
the  resolvent.  The  point  spectrum  of  B  consists  of  those  X 
for  which  (XI-B)  has  a  non-trivial  null-space. 

Theorem  1:  (i)  The  resolvent  of  Ap  is  {X  €  ff_Y:  det  A(X)  f  0} 

(ii)  The  point  spectrum  of  A^  consists  of 

{X  e  C_Y*.  det  A (X )  =  0),  and  if  det  A(-y)  *  0,  also 
X  =  -y. 

(iii)  The  spectrum  of  A^  contains  QNffy. 

(iv)  AD  has  compact  resolvent,  i.e.  (Ag-Xl)  1  is  compact  for 
all  X  for  which  it  is  a  bounded  linear  operator  on  C^. 
As  a  consequence,  the  generalized  eigenspace 
^(X)  Uj=1>K(AD-XI)k)  is  of  finite  dimension, 

denoting  the  null-space. 

For  (iv),  use  the  compactification  [-*,0]  *  ( -<■» , 0 )  u  {-°°},  the 
Arzela-Ascoli  Theorem,  and  a  well-known  theorem  on  projections 
(found,  for  example,  in  Hille  and  Phillips  [10,  p.  182]). 

Using  the  spectral  properties  of  A^  it  is  possible  to  obtain 


t 


6 


some,  but  not  all,  of  the  spectral  properties  of  Tp(t).  From 

Hille  and  Phillips  [10,  p.  467]  we  have  Pa  (TQ  (t)  )H  0}  =  exp (Pa(Ap) • t) , 

where  Pa  stands  for  the  point  spectrum.  For  the  rest  of  the 

spectrum  we  must  study  the  particular  equation  Dyt  =  0,  rather 

than  relying  only  on  abstract  results  for  semi-groups. 

Define  the  difference  operator  DQ :  ■*  ]Rn  by  Dp<t>  =  <J>(0), 

and  denote  by  Tn  (t)  the  semi -group  for  the  difference  equation 

0 

y(t)  =  DQyt  =  0  defined  on  the  Banach  space  €  C^:  D0<J>  =  0} 

with  the  |-|  norm.  As  an  aside,  note  that  |  Tn  ( t )  |  Y  <  e  1 4>  |  v 

r  uo  1  -  ' 

for  all  <J>  €  Define  the  projection  'Pq :  &->■  by 

=  MO  -  MO). 

Lemma  2:  Tp(t)  =  Tp  (tjH’p  +  U(t)  with  U(t)  completely  continuous 
on 

Let  now  |  •  |  also  stand  for  the  operator  norm  on  L(_#)  = 

(the  space  of  all  bounded  linear  operators  on  .  Using  Lemma  2 

and  arguments  involving  the  so-called  "essential  spectrum",  one  can 
prove,  as  was  done  by  Hale  [6,  p.  285]  and  Henry  [9,  p.  117]  for 
finite  delay: 

Theorem  3:  Let  aD  *  max [- Y, sup{Re  X:  det  A(X)  *  0}].  Then  for 
all  a  >  there  exists  K  =  K(a)  such  that  |Tp(t)|  <  Keat  for 

all  t  >  0.  One  calls  aQ  the  order  of  the  semi-group  Tp(t). 

Theorem  4:  For  any  -r  <  a  <  8  <  aD  the  set 

A  «  (X :  o  <  Re  X  <  8,  det  A(X)  ■  0)  has  only  finitely-many  points. 

For  Theorem  4  a  more  specific  reference  is  Hale  [6,  p.  309]. 


f 
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Using  knowledge  of  the  linear  problem  Dy^  =  0  we  can  discuss 
the  nonlinear  problem  (3)  by  obtaining  the  variation  of  constants 
formula  found  below  in  Theorem  5.  Let  X(t)  be  the  fundamental 
matrix,  i.e.  the  solution  in  ]RnXn  of  the  equation  DXt  =  I  for 

t  >  0  with  initial  data  XQ  given  by  XQ(0)  =  I,  XQ(0)  =  0  for 

9  <  0. 

Theorem  5:  The  general  solution  y^  €  of  the  inhomogeneous 
equation  Dy^  =  h(t),  y^  =  <P  €  Cy  with  h(t)  continuous  for 
t  >  0  and  with  D4>  =  h(0)  is  given  by 

ft  + 

(4)  xt  -  XQh(t)  =  TD(t)(<j>-X0h(0))  -  [dsTD(t-s)X0]h(s). 

Let  us  assume  now  that  G(t,0)  =  0,  D^GCt.O)  =  0,  and  that 

G(t,*)  depends  weakly  on  the  value  of  <f>  (0) ,  specifically  in  the 

sense  that 

(5)  G  C  t ,  4*  ±  XQb)  =  G(t,<t>),  for  all  <P  €  and  b  €  ]Rn. 

Define  a  new  space  PC^  =  ©  (the  span  of  the  columns  of  XQ) 

with  norm  |4  +  Xgb|.y  =  ) 4>  |  ^  +  |b|  n,  as  in  Hale  and  Martinez- 

1R 

Amores  [8].  Using  estimates  on  X  C  * )  and  the  measures  dsX('-s), 
as  in  Hale  [6,  p.  303],  the  space  PC^,  the  variation  of  constants 
formula,  and  Gronwall's  inequality,  one  can  justify  linearization 


using 


f 
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Theorem  6:  Assume  G  €  C*  0R+xCY  ;]Rn)  ,  G(t,0)  =  0,  D^G(t,0)  =  0, 
and  D^G(-,*)  is  bounded  on  all  sets  1R+  *  B,  B  bounded  in  C^, 
and  that  G  satisfies  (5).  Assume  that  |G(t,4>)|  <  Mg(|<|>|  ) 
where  M  is  a  positive  real  number  and  g  is  continuously 
differentiable  with  g(0)  =  0  =  g'(0).  If  aD  <  0  then  0  € 
is  uniformly  asymptotically  stable  for  equation  (3). 


Theorem  7  (Instability):  Make  the  same  assumptions  on  G  as  in 
Theorem  6.  If  a^  >  0  then  0  £  is  unstable  for  equation  (3). 

We  have  restricted  ourselves  to  D  of  the  form 

f0  Yg  ,0 

D$  =  *(0)  -  er  A(0)4>(6)de  with  |A(0)|d0  <  ». 

J  -00  j  -00 


All  of  the  above  results  can  be  achieved  when  one  allows  point 

- yr  rO  Q 

delays,  i.e.  for  =  <£(0)  -  ^^(-r^)  "  ®  A(9)^(9)d0 

J  -00 


fO 

with  IjJAjJ  +  I  |  A ( 0 )  |  d.Q  <  00 ,  as  long  as  one  assumes 


0  <  rl  <  r2  <  ’ ’ ’  (°r»  more  generally,  that  D  is  "atomic  at  zero"). 

We  discuss  the  case  with  point  delays  in  a  paper  currently  in 
preparation . 

Now  we  can. show  that  the  threshold  phenomenon  for  equation  (1) 
is  equivalent  to  the  question  of  the  stability  of  equation  (2). 

Since  equation  (2)  is  an  example  of  equation  (3)  for  which  G 
satisfies  the  assumptions  of  Theorem  6,  the  stability  of  (2)  can  be 
discussed  by  examining  the  stability  of  the  linear  difference 

r° 

operator  D(k)  given  by  D(k)$  *  (0 )  -  k  B(0)<)>(0)d0 ,  as  long 

f°  _Ye 

as  we  assume  that  e  |B(0)|d9  <  00  for  some  constant  Y  >  0. 


St 

I 
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/t  0  \ 

If  S  =  1/f  B(0)d0J  then  it  is  easy  to  see  that  (a)  k  >  S* 

'  -  00  ' 

jlj 

implies  a^^)  >  0>  and  (b)  k  <  S  implies  ap(k)  <  °*  Using 
f0  , 

S(t)  =  SQ  +  v(t+9)d9,  where  v(t)  =  ^  S(t)  and  v  solves  (2), 

along  with  Theorems  (6),  (7),  we  can  interpret  the  threshold 
phenomenon  as  a  particular  case  of  results  in  the  stability  theory 

of  functional  difference  equations:  For  k  <  S  ,  and  thus  ap(k)  < 

it  is  dear  that  S(°°)  =  +  0(|<f>|  ),  where  4>  (9 )  =  v(9)  =  ^  S(9) 

for  9(-°°,0].  For  k  >  S  ,  and  thus  aD(k)  >  t^le  semi-group 
Tp(t)  has  eigensolutions  e^*^,  where  <|>^(9)  =  e^9  is  in 

for  some  values  of  X  with  Re  X  >  0 .  In  fact,  X  =  is  such 
a  value!  This  latter  property  follows  from  the  fact  that  Tp(t)  is 
a  positive  operator  on  the  ordered  Banach  space  C  ,  since  A(9)  >  0 

■">  an 

9  €  ( -°» , 0 ] .  Let  <i>  =  <f>  ,  for  notational  convenience.  For  the 

epidemiological  problem,  the  initial  history  <t>  satisfies 
4>  (9 )  <  0,  9  £  (-“,0],  and  <f>  ( • )  f  0,  so  that  there  is  non-zero 

/v  /v 

projection  of  <)>  onto  the  subspace  [<t>]  of  spanned  by  <t> . 

From  this,  it  follows  that  lim^^S,,,  (<}> ; Sq)  <  Sq  whenever  Sq  >  S* , 

by  using  the  general  saddle-point  theory,  as  in  Hale  [4,  p.  157] 
or  Henry's  forthcoming  book,  Geometric  Theory  of  Partial  Differential 
Equations . 

Diekmann  [3]  has  allowed  spatial  effects  in  an  S  -*■  I  model 
to  arrive  at  the  equation  |-r  S(t,x)  =  S(t,x)*j  j  B  (9  ;  x  ,5  )S  (t+9  ,£  )  d£d9 

in  some  region  fi  c  ]Rm.  If  x  is  the  ordered  Banach  space  C(ft) 
we  can  re-write  this  model  as  a  functional  difference  equation  in 


the  ordered  Banach  space  Cy  ~  :  (-“,0]  Xj  <p  continuous  and 

Y0 

there  exists  11%.  <-6  ^  (0 ) )  •  If  we  restrict  the  linear  operator 
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D:  Cy  -*■  X  to  the  form  D<|>  =  4>(0) 


A  (6  )<{>  (9 )  d9  ,  with  the  map 


<J>  i— >- 


Y0A(9)<|>(0)d0  : 


X  being  completely  continuous  and 


order-preserving,  then  we  get  all  of  the  above  results,  especially 
Theorems  6  and  7.  From  this  we  can  get  threshold  results  for  the 
spatial  model  that  are  sharper  than  those  of  Diekmann  [3],  without 
the  monotonicity  assumption  of  Thieme  [14,  p.  103,  p.  94]. 

For  the  S  -*■  I  -*■  S  model,  where  infected*  can  recover,  it 
should  be  possible  to  prove  the  existence  of  periodic  solutions, 
via  Hopf  bifurcation,  when  the  spectrum  of  the  difference  operator 
D(k)  given  by 


D  (k)<J>  =  $(0) 


r 

J  -c 


eY0A(0)<j>(e)d6 


c 


eY0C(0)<t>  (0)d© 


depends  appropriately  on  the  parameter  k.  Smith  [13],  among  others, 
has  investigated  the  existence  of  periodic  solutions  above  the 
threshold. 
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